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ABSTRACT 

We point out a natural mechanism for quenching of star formation in early-type galaxies. It au- 
tomatically links the color of a galaxy with its morphology and does not require gas consumption, 
removal or termination of gas supply. Given that star formation takes place in gravitationally unstable 
gas disks, it can be quenched when a disk becomes stable against fragmentation to bound clumps. 
This can result from the growth of a stellar spheroid, for instance by mergers. We present the con- 
cept of morphological quenching (MQ) using standard disk instability analysis, and demonstrate its 
natural occurrence in a cosmological simulation using an efficient zoom-in technique. We show that 
the transition from a stellar disk to a spheroid can be sufhcient to stabilize the gas disk, quench star 
formation, and turn an early-type galaxy red and dead while gas accretion continues. The turbulence 
necessary for disk stability can be stirred up by sheared perturbations within the disk in the absence 
of bound star-forming clumps. While other quenching mechanisms, such as gas stripping, AGN feed- 
back, virial shock heating, and gravitational heating are limited to massive halos, the MQ can explain 
the appearance of red early-type galaxies also in halos less massive than ^ IO^^Mq. The dense gas 
disks observed in some of today's red ellipticals may be the relics of this mechanism, whereas red 
galaxies with quenched gas disks could be more frequent at high redshift. 
Subject headings: galaxies: elliptical and lenticular, galaxies: evolution, galaxies: formation 
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1. INTRODUCTION 

In the low redshift Universe, the variety of mor- 
phological t ypes of galaxies is accomp anied by a color 
bimodality ("Rob erts fc Havnes I 119941 ). Spiral galax- 
ies have low bulge-to-disk ratios and are actively star 
forming, with blue colors. Early-type galaxies (ETGs), 
which are spheroid-dominated lenticulars and ellipticals, 
are red, non- or almost non-star-forming. The for- 
mation of disk galaxies could be explained by a va- 
riety of processe s, from smooth gas accretion and in- 
ternal evolution (lEl megreen ct al. 2005; Bournaud et al.l 
I2007at iGenzel et al, 2008: Dekel et al. 2009a.b.l to the 
early merging of gas-rich galaxies (iRobertson et all 
l2006l:lSpringel fc Hernauistll2Cl05tlGovernato et al.ll2007D . 
There is a larger consensus that the transformation 
of spiral disks into elliptical galaxies is largely driven 
by major and minor galaxy mergers dBarnesi il992 : 
INaab fc Burkertl [2001 feournaud et^ 120051 . I2"007bl ). 
The morphological transformation of spirals into ellip- 
ticals thus has to be accompanied by a "quenching" of 
their star formation, which shifts them to the Red Se- 
quence. 

This quenching of star formation can be obtained by 
consuming the cold gas reservoir. However, simula- 
tions suggest that this might not be the c ase for typ- 
ical mergers ()di Matteo et al. 1 120071 I2008D . In stan- 
dard pictures, some extra quenching processes are then 
needed to completely remove the cold gas reservoir 
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of these post-merger galaxies and to shut down any 
fresh supply of cold gas. Such processes include (i) a 
cluster environment in which the cold gas is stripped 
and no cold gas can be found in the neighborhood 
(but see iTonnesen fc BrvanI 120091 : iSalome et all 120081 : 
'Revaz et al.ll2008l ). (ii) heating by virial shocks of gas in- 
falling in a massive dark raatter halo (Birnboim & DeM 
2003: iDekel fc Birnboi mI (20061: iBirnboim et al. I 120071 
j Cattaneo et al.l .2008 : Ocvirk et al. I 120081: iKeres et all 
I2009D . fiii)" energy feedback from an Active Galac- 
tic Nucleus (AGN}, which could itself be trigger ed by 
galaxy mergers (jHoDkins ct al. 2005; Di Matte o et al 



20081 iCatta neo et al. 2005, 2006. 20091 : 1 Johansson et al 
2009al: iScha win ski et al., 2QM) (iv) and gravita- 
tional lieating__o^_a^_gas disk by clumpy mass in- 
fall (IDekel fc B irnboimI 120081) and minor mergers 
pohansson et al . 2009b). While AGN feedback could 
participate in the reddening of massive galaxies and 
nuclear activity i s frequently observed at high redshift 
(jNesvadba et al. 1 [2008) , there is a lack of direct evi- 
dence that this is the looked-for quenching mechanism 
- AGN jets could be too coUimated to affect most cold 
gas disks, and AGNs could theoretically trigger the star 
formation as well. Gas heating by AGN feedback and 
virial shock is furthermore lim ited to halos above ^10^^ 
Mq (Dekel fc Birnboimll2008h . 

The mechanisms listed above therefore struggle to 
explain the reddening of fie ld ETGs in halos belo w 
this critical virial in ass (e.g.. iWeinmann et al.l l2006al lbl: 
iHopkins et ani2008f ). In this paper, we propose a new 
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mechanism called "morphological quenching" (MQ), in 
which star formation is quenched even for ETGs in halos 
of moderate masses. Contrary to the mechanisms pro- 
posed so far, we show that an ETG can spontaneously 
turn red without having its cold gas removed, and could 
even contain amounts of cold gas similar to those of some 
spiral galaxies. This gas would be stabilized against star 
formation by the morphological transition from a rotat- 
ing stellar disk to a pressure-dominated spheroid, which 
induces a steeper potential well and reduces the disk self- 
gravity. Observations of^disk galaxies with various bulge 
fractions pCennicutt Ill989i ) had already shown that bulge 
dominated galaxies form stars less efficiently than disk- 
dominated ones. We propose that this morphological 
quenching of star formation could participate in the red- 
dening of elliptical and lenticular galaxies. 

This mechanism is found to naturally occur in numer- 
ical simulations of galaxy evolution in the A-CDM con- 
text. We use a new technique for cosmological zoom-in 
simulations, in which the cosmological history of a given 
galaxy (merger tree and smooth accretion) is recorded, 
and re-simulated at high-resolution. Our simulations 
suggest that weak gravitational instabilities in the gas 
disk (spiral arms) can trigger the turbulence necessary 
for stability while not allowing dense bound clouds to 
form, preventing efficient star formation. 

Morphological quenching is in agreement with obser- 
vations showing that some Local red ETGs host sig- 
nificant amounts of atomic or molecular gas with high 
densities, sometimes in amounts similar to those of 
the Milky- Way gas disk, but with only very inefficient 
to no star formation at al l and red integrated colors 
jMorganti et al. 2006; Croc ker et all l2009l:T^ung et all 
[2009bt iDonovan et al.ii2009l) . This process could be even 
more important at hig h redshift where a red sequence is 
likely already in place (jBell et al. Il2004t ICirasuolo et al. I 
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I2007D while gas fractions in early-type galaxies are ex- 
pected to be higher since gas infall rates are high and 
most clusters are not virialized yet. 

Our simulations also indicate that a reddened ellip- 
tical galaxy could turn back again into a blue star- 
forming disk galaxy provided that it re-accretes sufficient 
amounts of cold gas. This transition back to the blue se- 
quence requires high rates of cold gas accretion, probably 
making this process infrequent at z < 1. 

The MQ mechanism is described in Section 2. Cosmo- 
logical models illustrating the formation of ETGs with 
morphological quenching are presented in Section 3, and 
the stability of their gas disk against star formation is 
studied in Section 4. Implications for various observa- 
tions are discussed in Section 5. 

2. MORPHOLOGICAL QUENCHING 

Gravitational instability is a necessary con dition for ef- 
ficien t star formation i n a galactic disk (e.g. , lElmegreen I 
[200l iLi et al. I [200l iKawata et aT~l [20071 ) The insta- 
bility of a single-component, rotating disk at radius r 
is characterized by the Toomre Q parameter (jToomre I 
[19641 : 



(1) 



Here a is the one-dimensional velocity dispersion and S 
is the surface density of the disk, k is the epicyclic fre- 



quency, which is linked to the depth of the gravitational 
potential, and is of the order of the angular velocity fi. 
The factor e is 3.36 for a thin stellar disk and tt for a 
gas disk. Axisymmetric instabilities, namely rings that 
break into clumps, can grow in the disk if Q < 1. 

The instability of a two-component (gas and stars) 
disk is characterized by an effective To omre parameter 
(|Jog fc Sob mon 1 [l98l lElmegreenl[l99l) defined by 



as 
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where the Q parameters for the gas (g) and for the stars 
(s) are computed separately for each component using 
Equation [T] The coefficients are 
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gs — cTg/cTs and the dimensionless wave number 



where a. 

q = kus/ n is the one minimizing Q. These coefficients 
are of order unity in spiral galaxies. 

The stellar surface density Ss in Equation [1] refers to 
the stars that co-rotate in a relatively thin disk and not 
to the stars that populate the spheroids, bulge or halo. 
In particular, for an elliptical galaxy with all the stars 
in a pressure-supported spheroid, Qs oo . A lentic- 
ular galaxy with a thick, kinematically hot stellar disk 
embedded in a massive bulge has Qs » 1. 

The stability of a disk for clump-forming instabilities 
and subsequent star formation has indeed been shown 
to be properly characterized by the effective Q parame- 
ter both in observations (Wong & Blitz 2002, Leroy et 
al. 2008) and in simulations (Li et al. 2005). A disk 
with a sufficiently low gas surface density and/or a suf- 
ficiently high epicyclic frequency has Qg > 1, which, by 
itself, would have indicated stability for axisymmetric 
modes. However, for a given Eg and k, the presence of a 
stellar disk component enhances the self-gravity within 
the perturbation and can destabilize the disk, pulling 
the effective Q parameter to below unity. For exam- 
ple, typical low- redshift disk galaxies have CTs >> Cg 
as well as Eg >> Eg, so Qs could be comparable to 
Qg and both can thus play a similar role in determin- 
ing the disk instability. The same likely holds at high 
redshift, where Eg Es but fjg ~ as (see for instance 
iForster Schreiber et al.l l2009). The Q < 1 condition en- 
sure the formation of gas clumps in which new stars form, 
without implying that the pre-existing stars collapse in 
these gas clumps. 

Furthermore, a necessary requirement for star forma- 
tion is that the internal gravity in gas clouds overcomes 
the disruptive galactic tidal field. This requires a three- 
dimensional disk de nsity above the critical tidal density 
([Hunter et al.l[200l 
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Note that the relevant disk density involves both the 
gaseous and stellar components of the disk; the disk stars 
may help the self-gravity of the clumps overcome the ex- 
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ternal tidal forces even in cases where the disk gas is not 
dense enough for the purpose. 

The basic idea of morphological quenching is that the 
star formation in a gaseous disk could be severely sup- 
pressed even if the disk is relatively massive once it is 
embedded in an early-type galaxy, that is in an ellipti- 
cal or a bulge-dominated lenticular. A similar gas disk 
could be effectively star-forming in a spiral galaxy where 
most of the stellar mass is in a rotating disk. As a result, 
the former type of galaxies would naturally have red col- 
ors, while the latter would be bluer. The difference in 
disk stability arises from two main effects: (i) the higher 
concentration of the stellar mass in an ETG, affecting Q 
via K, and (ii) the contribution of the stellar disk to the 
self-gravity in a spiral, affecting Q via Eg. From another 
perspective, when the stellar disk becomes a spheroid in 
an ETG, the self-gravity of the gas disk that is left be- 
hind can fail to balance the disruptive tidal forces, which 
prevents the assembly of star-forming clumps. 

Bulge-dominated galaxies are therefore likely to form 
stars less efficiently than disk-dominated ones, as ob- 
served bv lKennicutt I ()1989f ) for disk galaxies with various 
bulge fractions. We here propose that this effect could 
actually lead to quenching of star formation in E and SO 
galaxies and make them red. Furthermore, in the follow- 
ing, we show in numerical simulations that morphological 
quenching can naturally occur in the cosmological con- 
text, and that gas disks with masses up to several 10^ Mq 
can be stabilized in ETGs. Our simulations suggest that 
a rather high turbulent speed ((Tg) can be maintained in 
a stabilized gas disk in an ETG (see Section 4), which 
provides the pressure support necessary for stability. 

3. NUMERICAL SIMULATIONS 
3.1. Method 

To study galaxy evolution in the cosmological con- 
text, the standard technique consists in running full cos- 
mological simulations with spatial and mass resolutions 
increased in a certain "zoom - in" area (|Katz fc W hite 
1991 iGove rnato ct al. 2004; Scmcl in fcCombed (20051 : 
Naab et al. . ,2007; Joung et a l. .200^ A constraint is 
that all particles ending-up in the target galaxy should 
be low-mass, high-resolution ones. This makes uneasy, 
if possible at all, to simulate galaxies with high accre- 
tion and merger rates at all redshifts, while more iso- 
lated galaxies or massive galaxies with early assemblies 
are more easily modeled this way. 

We here use an original approach where we extract 
merger and accretion histories (and geometry) for a given 
halo/galaxy in a A-CDM cosmological simulation (see 
Appendix A.l), and then re-simulate theses histories at 
much higher resolution, replacing each halo by a realis- 
tic ga containing gas, stars and dark matter (Appendix 
A. 2). The re-simulation is performed with a particle- 
mesh code where gas dynamics is modeled with a sticky 
particle scheme (Appendices A. 2 and B). This method 
takes into account not only every merger prescribed by 
the cosmological simulation, but also the smooth gas and 
dark matter infall along cold flows or cosmic filaments. 
A similar approach has been independently proposed by 
Kazantzidis et al. (2008), Read et al. (2008) and Villalo- 
bos & Helmi (2008), but these authors account only for 
the main mergers and not for the diffuse accretion, and 



only performed coUisionless simulations without any gas 
component, and sometimes even no baryonic component 
in satellite galaxies. Here we model both the main stud- 
ied galaxy and all companions with dark matter, stars, 
and gas particles, and include all the accreted mass, from 
major and minor mergers to smooth accretion (Appen- 
dices A. 2.1 and A. 2.2). This allows us to study galaxy 
evolution in a realistic cosmological context from a high 
initial redshift (z = 2 in the present paper) down to z = 
and account for the internal processes at high resolution: 
the maximal spatial resolution (and softening length) is 
130 pc, while the mass resolution is 1.4x10^ for stel- 
lar particles initially present in the galaxies, 2.1x10^ 
for gas particles and stellar particles formed during the 
simulation, and 4.4x10^ M0 for dark matter particles. 
The technique is described in detail in Appendix A. Its 
main advantage is a reasonable computer cost for all kind 
of galaxies, even in dense environments and with high in- 
fall rates at all redshifts. Large samples of simulations 
with this technique will be studied in future papers. 

In the present paper, we focus our analysis on the color 
evolution of a galaxy similar to the Milky Way in terms 
of mass at z = and which does not belong to a rich 
group or a cluster of galaxies. The host halo is shown 
with its cosmological environment at z = 2 and z = 
on Figure [H During the simulation, the halo mass grows 
from 2x10" Mg at z = 2 to 1.4x10^2 Mg at z = 0. This 
galaxy undergoes a first series of several minor mergers 
soon after z = 2, four of which have mass ratios between 
4:1 and 10:1. Then comes a quieter phase that is dom- 
inated by smooth dark matter and gas accretion, with 
an average gas accretion rate of 9 Mg yr~^, lasting from 
z ~ 1.2 down to z ~ 0.2. Finally, a major merger takes 
place at z = 0.2 with a mass ratio of 1.5:1. It is shortly 
preceded by a peak in the diffuse gas accretion rate (up 
to 70 Mg yr-i). 

A case without the final major merger but the same 
history for other mergers and diffuse mass accretion is 
also studied. In this modified simulation, the whole mass 
corresponding to the incoming galaxy (both dark matter 
and baryons) is removed. 

3.2. Results 

The morphological evolution from redshift 2 to red- 
shift is shown on Figure [2l The evolution of stellar 
and gaseous masses, star formation rate (SFR) and inte- 
grated B — V color are shown on Figures [3] and ID These 
quantities are integrated in a sphere of radius 25 kpc 
around the main galaxy (see Appendix A. 3 for details on 
the analysis technique). 

3.2.1. Early mergers and spheroid buildup 

The simulation starts with a gas-rich disk galaxy (see 
Appendix A. 2.1), forming stars actively in giant clumps, 
consistently with obser vations of high-red shift galaxies 
(jElmegreen et al.l 120071 : IGenzel et all '20081. A series of 
minor mergers takes place between redshifts 2 and 1, four 
of which have mass ratios between 4:1 and 10:1. These 
minor mergers transform the initial disk galaxy into an 
elliptical-like one, i.e., a st ellar spheroid with a low rota- 
tional support {V/a <1) (jBournaud et al.l l2007bl ) . The 
most important mergers occur at z ~ 1.5, and by z ~ 1 
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Fig. 1. — Large-scale dark-matter distribution in the cosmological simulation. Left: at z = 2 in a box of side 2 Mpc. Right: at 2 = in a 
box of side 4 Mpc. The halo on which we zoom in is marked by an arrow in the left panel and is at the center in the right panel. The 2 = 2 
snapshot shows neighboring halos that will merge with our halo at z ~ 1.5. The z = snapshot follows a major merger that occurred at 
z = 0.2. The halo is continuously fed by gas streams that follow the dark-matter filaments seen in these maps. 
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Fig. 2. — Evolution of the distribution of gas (top panels) and stars (bottom panels). The snapshot at t = 13.5 Gyr is shown both 
for the original simulation and for its variant performed without the final major merger. The box size is 40 kpc in physical units. The 
sequence of snapshots highlights the three phases of evolution. Until t = 6.5 Gyr, there is a gradual build-up of a stellar spheroid due to a 
suite of mergers. At t = 10.5 Gyr, the system is toward the end of the MQ phase, showing a large gas disk and a massive stellar spheroid. 
At t = 13.5 Gyr, after the major merger, the galaxy is an elliptical, but if this merger is eliminated, it becomes an ETG with massive disk 
and spheroid components. 



the galaxy has an elliptical shape, even if still interacting 
with a few lower-mass companions. 

This morphological transformation goes along with an 
episode of active star formation, reaching 50 Mq yr^^ 
(right panel of Figure [3]). At the same time, external gas 
is infalling along a filament, which maintains a high gas 
fraction (the right panel of Figure [3|) and sustains a high 
star formation rate during 2-3 Gyr. 

As for the integrated colors, their evolution during the 
first 1.5 Gyr is sensitive to the assumptions made on the 
star formation history of each galaxy before its entry into 
the virial radius (see Appendix A), which can be seen 
as an initial scatter in Figure [H Under most realistic 
assumptions, the initial galaxy is blue {B ~ V < 0.5). 



The starbursts taking place during the minor mergers 
rapidly wash out any initial difference, and result in a 
very blue galaxy with B — V = 0.28 a.t t — 6 Gyr. 

3.2.2. Morphological quenching phase 

After the initial mergers, a much quieter period lasts 
from z 2± 1 to z = 0.2, during which no merger takes 
place other than with tiny satellites. This phase is qui- 
escent in terms of merger events; however, diffuse gas 
is continuously infalling at a steady rate; the average 
gas accretion rate across the virial radius of the studied 
galaxy during this phase is 9 Mq yr~^. 

The infalling gas adds to the gas mass left over after 
the initial mergers, so that the cold and dense (p > 0.01 
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Fig. 3. — Evolution of the stellar and gas components in the central sphere of radius 25 kpc. Left: mass in stars (solid, blue) and in 
baryons (dashed, red). Right: the fraction of gas with density above 0.01 Mqpc"^ with respect to the total baryonic mass and the star 
formation rate. The shaded area highlights the morphological-quenching phase, where the SFR is suppressed and the stellar content is 
almost constant while the gas content is increasing gradually due to continuous accretion. The earlier and later phases are characterized 
by high gas fraction and high SFR. The growth in the early phase is dominated by a sequence of minor mergers, and in the late phase by 
a major merger. 
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Fig. 4. — Evolution of the B — V color, showing the three phases 
as in Figure 3. The solid (dark grey) curve corresponds to the 
fiducial simulation. The galaxy is blue during the early and late 
phases, and is reddened by half a magnitude during the MQ phase. 
The dark shaded area along the curve represents the range of colors 
obtained for different star formation histories (Appendix |XJ . The 
dash-dotted (blue) line corresponds to a star formation threshold at 
0.25 M0pc~'' (Section 4.3), showing a somewhat stronger quench- 
ing toward the end of the MQ phase. The dashed line (red) refers 
to the case where the final major merger is eliminated, showing 
only little differences. The dotted curve shows the evolution after 
2 = 0, where the system is about to enter a new quenched phase. 
The reddening associated with MQ is a robust phenomenon. 

Mopc"'^) gas mass in the studied galaxy varies from 
5.1x10^ M0 at < = 8 Gyr to 1.3 x 10^° Mq at t = 11 Gyr 
(the left panel of Figure E]). During the quiescent phase, 
the galaxy thus remains relatively gas-rich, with a dense 
gas fraction increasing from 4% to 8% (with respect to 
the total baryonic mass inside a 25 kpc radius, see Fig- 
ure [3]). 

This gas settles into a disk that slowly grows with time 
and acquires a tilt of almost 90 degrees with respect to 
the initial disk plane (Figure [5] and O because the main 
filament providing the infalling gas is almost perpendicu- 
lar to the initial disk plane. This gas disk grows from the 



Fig. 5. — The galaxy during the morphological quenching phase 
at t = 10 Gyr. The box side is 30 kpc. The stellar mass distribution 
is represented in yellow-red colors and the gas distribution in white- 
purple colors. The gas disk undergoing morphological quenching of 
star formation can be seen inside the spheroidal stellar component. 

inside out: the density of the inner regions (r < 2.5 kpc) 
does not evolve much during the quiescent phase while 
the accreting gas settles in the outskirts of the disk . 

During this period, star formation is not completely 
suppressed but the star formation rate remains nearly 
constant at around 1.5 Mq yr~^, which is a low level 
considering the gas content of the galaxy, as shown in 
Figure [3l between z = 0.6 and z = 0.2 the gas fraction 
increases by a factor of 2 while the SFR remains low. 
For such gas fractions, and more importantly such gas 
surface densities (see also Section 4.1 and Tabled]), star 
formation rates above 10 Mq yr~^ could have been ex- 
pected: for instance the galaxy at the beginning of the 
MQ phase has the same gas mass and gas surface density 
as the disk galaxy at t = 3.5 Gyr but the latter forms 
stars at a rate of 9 Mq yr"-"^ (this will be discussed in 
more details in Section 4.1). 
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The low SFR during the MQ phase aUows the forma- 
tion a low-mass young stellar disk that never exceeds a 
few percent of the total stellar mass: on average 2%, and 
at most 6% at t=ll Gyr which is the end of the quiescent 
phase. Thus, most of the stars of this galaxy still belong 
to the spheroid formed by the series of minor mergers at 
I < z < 2. The galaxy therefore remains elliptical-like, 
with only a low-mass disk component, as in many real 
elliptical galaxies (see lKrainovic et al. II2008I ). 

As for the color evolution, the galaxy rapidly reddens 
after the z >1 mergers that converted it into an elliptical- 
like system (FigureH]). During more than 4 Gyr, it main- 
tains a red color with B — V = 0.7-0.8 most of the time, 
w hich is consistent with colo rs of local ETGs according 
to [Roberts fc Havnesi (jl994[ ). The integrated color is 
indeed dominated by the contribution of the relatively 
old stars formed before/during the series of mergers. Be- 
cause of the low star formation rate, young stars with 
bluer colors are not numerous enough to impact the in- 
tegrated color, even though their contribution is max- 
imized by neglecting dust extinction in our color esti- 
mates. 

The same color evolution is also found in other color 
indices such asU ~ B and U — V (see Appendix C). 

This elliptical galaxy remains constantly red during the 
quiescent phase in spite of a significant gas fraction left 
over after the initial mergers and constant infall of gas at 
a steady rate. Indeed, no mechanism removing the cold 
gas and/or terminating the accretion of fresh gas (like 
AGN feedback or virial shocks) occurs in our models. As 
a consequence, this red elliptical hosts a dense and mas- 
sive disk of cold gas (Figure [5]): in average 5 x 10^ of 
cold gas with typical surface densities of 5-10 M0 pc~^. 

This simulation shows that the effect of morphological 
quenching can be strong enough to turn an ETG red 
without requiring the removal of cold gas reservoirs and 
the termination of gas supply. 

We will show in Section 4 how the low efficiency of star 
formation in this gas-rich red elliptical is connected to the 
presence of a stellar spheroid instead of a disk-dominated 
galaxy. We will also show that the gas disk is stable 
against star-forming instabilities, but unstable against 
non-axisymmetric modes which trigger the turbulence at 
the level required for Q > I stability. 

3.2.3. Late blue phase 

The studied galaxy finally undergoes a major merger 
between z = 0.2 and z = 0, together with an increase of 
the diffuse gas accretion rate. 

The combination of the major merger and the in- 
creased gas content allows a large starburst to take place, 
with a peak SFR of up to 35 M© yr~^. As a consequence, 
the galaxy rapidly evolves toward blue colors. 

However, this blue color at z = is only a transient 
phase related to the merger-driven starburst. We have 
let this galaxy evolve after z = 0, and find that the star 
formation rate steadily decreases from 30 Mq yr~^ at z = 
to 10 M0 yr"-'^ after 2 Gyr. The galaxy reddens again, 
reaching B — V > 0.5 within 2 Gyr (Figure [3|). This red 
early-type galaxy still contains a significant fraction of 
cold gas (6% of the baryons), both because some gas is 
left over after the final merger, and some infalling gas 
that entered the virial radius at low redshift continues 



to feed the system (even if we neglect gas accretion at 
the virial radius after 2 = 0). We thus find another 
ETG turning red spontaneously, under the effect of the 
MQ mechanism and without any mechanism removing 
its cold gas reservoir and/or stopping the infall of fresh 
gas. 

To better understand the relation between morphol- 
ogy and color evolution, we have performed the same 
simulation without the final z ~ 0.2 major merger - ear- 
lier mergers and diffuse accretion remaining unchanged. 
While cosmologically realistic since directly motivated by 
a A-CDM simulation, the fiducial case may be less likely 
than this version without the final major merger, since 
the majority of massiv e galaxies do not undergo a major 
merg er at z < 1 (e.g.. iGenel et al.l 120081 : iStewart et al. I 

Hooi. 

With this model, we find that it is the accumulation 
of cold gas accretion more than the final major merger 
itself that causes the return to blue colors at z ~ 0.2 (Fig- 
ure |4]). Indeed, even without the final major merger, the 
galaxy undergoes a phase of intense star formation at low 
redshift, with a star formation rate reaching 25 M0 yr"-'^ 
at z = 0. The most important difference induced by re- 
moving the final merger lies in the morphology of the 
galaxy at z = (Figure [5]). With the final major merger, 
the galaxy is elliptical, and we find that its blue color 
is only transient and rapidly turns red under a MQ af- 
fect (see above). Without the final merger, the galaxy is 
more disk dominated, forming stars more durably; it de- 
velops spiral arms that host efRcient star formation and 
are visible in the z = stellar disk (Figure^ see also Sec- 
tion 4). Hence, the gas disk that was stabilized against 
star formation during the quiescent red phase can be- 
come actively star-forming again if the accumulated gas 
mass increases too much. 

3.2.4. Summary 

Our main findings from the simulation so far can be 
summarized as follows. First, a spheroid-dominated 
galaxy naturally becomes red and dead despite the pres- 
ence of a relatively massive gas disk. There is no need 
for an external mechanism that removes the present gas 
or terminates the supply of fresh gas. Second, a red 
spheroidal galaxy can turn back into a blue system com- 
posed of a bulge and a disk. In our current simulation, 
this happened as a result of a continuous growth of the 
disk over several Gyr through smooth gas accretion. The 
rather smooth gas streaming into a galaxy could thus sig- 
nificantly affect the morphology and the color of a galaxy. 
This is in addition to the effects of mergers, which tend 
to build spheroids and thus quench instability and star 
formation in gas disks. 

4. DISK STABILITY 

4.1. Stability analysis 

We find in our simulations that a dense gas disk can be 
stable against star formation when it is embedded in an 
early-type galaxy instead of a stellar disk. We will show 
in the following that a low gas surface density is not the 
reason for stability in the red ETG, and that the same 
disk would form stars efficiently if it was embedded in a 
stellar disk instead of a stellar spheroid. 

The cold gas disks in the red elliptical galaxy at t = 
9 Gyr and in a blue spiral at t = 3.5 Gyr have similar 
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TABLE 1 

Comparison of galaxy properties at different instants in the simulation. The values corresponding to 13.7 
Gyr are taken from the simulation without the final major merger, so that a gas disk can be identified. 



SFR(MQyr-i) B-V Mg (lO'' Mq) o-g (km s"!) cts (km s"!) < Qg > < Qs > < Qtot > 

t = 3.5 Gyr 9.1 0.39 3.96 25 43 8.7 4.6 2.3 

t = 9 Gyr 1.5 0.74 4.03 16 26 10.9 10.1 5.2 

t = 10 Gyr 1.4 0.76 6.14 15 19 9.1 7.6 4.2 

t = 13.7 Gyr 27.2 0.30 11.23 22 33 8.2 2.3 2.2 




Radius (kpc) 



Fig. 6. — Surface-density profiles of cold disk gas in the early 
blue phase (solid, blue) and in the MQ phase (dashed, red). The 
times are t = 3.5 and 9 Gyr respectively. The two disk profiles are 
similar, while the corresponding colors and SFRs are very different. 

masses, radii and densities (Figure [S]), but the former 
forms stars at a low rate, while the latter forms stars 
efficiently ^. The total (gas+stars) Toomre parameter 
for these two disks has an average value of 5.2 in the 
former case (red ETG), versus 2.3 in the latter one (blue 
disk galaxy) (Table [T|) . The difference in the stability of 
these gas disks and their efficiency of star formation does 
not result from different gas masses or densities, but from 
a different depth of the gravitational potential measured 
by K ^ and from a different contribution of stars Sg to 
the total disk density. 

Note that the disk instability and the high level of 
star formation at t = 3.5 Gyr are not artifacts due 
to relaxation of initial conditions. Indeed, the galaxy 
has been evolved in isolation for 500 Myr before being 
introduced in the simulation, which is a long enough 
time for the initial disk to acquire a realistic struc- 
ture in the global galactic potential (see more details in 
iMartig fc Bournaud ,.200a ). 

Maps of the Toomre parameter for the gas component 
(Qg), the stellar component (Qg) and the combined pa- 
rameter Q are shown in Figure [7] for the blue galaxy 
at t=3.5 Gyr, for the inefficiently star-forming red ETG 
at <=10 Gyr (i.e. during the MQ phase), and the ac- 

^ If star formation is calibrated so that the spiral galaxy follows 
the Schmidt-Kcnnicutt relation, the galaxy during the MQ phase 
then lies below this relation, but within the observed scatter. In- 
deed, this scatter is rather high, and allows SFRs to differ by a 
factor of 10 or more for a given gas surface density. On the other 
hand, changing the SFR by such a factor is enough for a galaxy to 
turn red. 

K is computed in radial bins of 100 pc from a smoothed version 
of the rotation curve. 



tively star- forming disk galaxy at t—13.7 Gyr in the case 
without the final merger so that a disk can be identified. 
These maps confirm that the star-forming disks have a 
low Q parameter, with averaged values of Q around 2 
in the star-forming regions, and local overdensities with 
Q < I in clumps and spiral arms. Regions with Q ~ 1 are 
unstable against axisymmetric instabilities, reach high 
gas densities, and form stars actively. 

The gas disk in the red ETG differs by having some- 
what larger values of Qg in its dense regions, because 
of the change in k, and more importantly by having 
Qs >> 1, because of the lack of stellar contribution Eg 
to the disk density and self-gravity. As a result, the com- 
bined Q parameter remains significantly above 1 every- 
where in the disk, preventing axisymmetric instabilities, 
even though non-axisymmetric armlets are present in the 
disk (Figure [HI). During the MQ phase, the gas thus does 
not reach high densities, and the star formation efficiency 
is low: considering the global cold gas mass present in 
the galaxy, the SFR could have been expected to be ten 
times higher. 

During the red phase of our simulation, the gas mass 
and surface density increase by a factor of 3. Yet the 
SFR remains constantly low and the galaxy remains red: 
the stability parameter Q indeed remains large, and no 
dense gas clumps form. Only when the mass of the cold 
gas disk has increased above 9x10^ Mq , the disk becomes 
unstable again, so that the star formation increases and 
the galaxy turns back blue: we find an average Q ~ 3 
with the densest regions reaching Q = 1 at t = 11 Gyr, 
when the galaxy begins to turn back to blue colors. This 
is illustrated by the simulation where the system evolves 
without the final major merger: the accreting cold gas 
continues to settle into the disk, increasing its surface 
density. Grand-design spiral arms develop (Figure [5]), 
and efficient star formation occurs in these arms. This 
results in a blue galaxy with a well-defined spiral struc- 
ture. Here again we note that the star formation activity 
and the morphological evolution remain connected: when 
the galaxy turns back blue, it forms a massive stellar disk 
and rapidly takes a spiral disk-|-bulge morphology. 

4.2. Turbulence in the stable phase 

We have shown above that an early-type galaxy tends 
to stabilize its gas disk against star formation. However, 
this analysis holds only if the gas turbulence provides the 
necessary pressure support. If the velocity dispersion a 
was reduced in early- type galaxies, their gas disks could 
be Jeans-unstable (low Qg) in spite of the absence of a 
massive stellar disk. The velocity dispersion of the gas 
disk particles in our simulation remains about similar 
between the star-forming and non-star forming phases, 
where it lies on average around 15 km s~^ (Table [Ij. 
For star-forming galaxies, Tamburro et al. (2009) indi- 
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f = 3.5Gyr f = 10Gyr f=13.7Gyr 

Star-forming clumpy Quenched gas disk Star-forming 
high-z disk in a red ETG low-z spiral 




Fig. 7. — Gas surface density and Q parameter in the disk. Shown are gas surface density (top), Qg, Qs and the combined effective 
Q in the three phases of evolution. The box side is 30 kpc. During the MQ phase, the effective Q is above 2 across most of the disk, 
but below 3-4 in extended areas. This explains the stability of axisymmetric modes and the lack of bound clumps in this phase, while 
non-axisymmetric perturbations are clearly visible. In the early and late phases, Q drops below unity in several extended areas, consistent 
with the appearance of strong axisymmetric perturbations and bound clumps in these phases. The effective Q is similar to Qg in the MQ 
phase, where is rather high. In the unstable phases, especially the late one, the stellar disk component helps destabilizing the disk. 




Fig. 8. — Surface density maps of the gas disk during the MQ phase. The snapshots are 33 Myr apart t = 10 Gyr. The box side is 15 kpc. 
The spiral armlets are strongly sheared and short-lived. For instance the dense structure A is just starting to form on the first snapshot, 
and already dissolving on the last one; structures B and C are not seen yet on the first snapshot, and already weaken on the last one. 
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cate a characteristic HI linewidth value of 8-12 km s~ 
in 11 spiral galaxies. This value includes turbulent ve- 
locity dispersion as well as thermal broadening, but tur- 
bulence is expected to dominate at least in the central 
(star- forming) regions of these galaxies. The velocity dis- 
persion of our sticky particle medium is slightly higher^ 
but includes all gas phases, so is consistent with obser- 
vations on the blue sequence. A more important wonder 
is the level of turbulence in early-type galaxies with low 
star formation efficiencies. 

The ISM turbulence in star-forming spiral galaxies is 
mostly triggered by gravitational instabilit i es and star- 
form a tion feedback (lElmeereen et al.ll200£ : iWada et al.l 
I2OOI iTasker fc BrvanI l200a lAgertz et al.l 120091 ). An 
equally high level of turbulence in an inefficiently star- 
forming disk could then seem unexpected. Several pro- 
cesses can nevertheless maintain the turbulence in the 
absence of star-forming instabilities: these could be (i) 
gravitational instabilities without sta r formation, (ii) ro - 
tational shear powering turbulence (j Shore "eraLll2n03h . 
(iii) interactions with small satellit es, dark subhalos , 
or more generally clumpy infall (Dcke l et al.l l2009bf ). 
and (i v) the magnetorotati onal instability in magnetized 
disks (jSellwood k Balbus"! [I999l) . Mechanism (ii) does 
not occur in our simulations since the critical Reynolds 
number above which rotational shear inst abilities might 
be able to grow is largely beyond reach (jJi et al.l 120061 : 
iLesur fcT^o ngaretti 200^. Mechanism (iv) is not in- 
cluded in our models, and mechanism (iii) cannot domi- 
nate as the mass infall on the studied red galaxy is quite 
smooth. Mechanisms (i) should then be the main source 
of turbulence. 

Gravity-driven turbulence does not necessarily require 
axisymmetric instabilities like dense clumps, and can ac- 
tually be triggered with Q > 1. A study of flocculent 
spirals bv lElmegreen et "al] (|2003f ) shows that turbulent 
motions are generated by non-axisymmetric instabilities 
(spiral arms). This turbulence then acts at smaller scales 
where it form overdensities (clouds), which in turn col- 
lapse if they exceed the critical density for Q=l. In 
this picture, the turbulent energy is mostly pumped 
by large-scale, non axisymmetric instabilities like spiral 
arms. These structures can form with values of Q much 
higher than 1 and trigger the turbulence, without involv- 
ing the collapse of gas clumps where stars would form effi- 
ciently. Density maps of the non-star-forming gas disk at 
t = 10 Gyr support this picture (Figure [5]): small spiral 
arms are ubiquitous. They have small pitch angle typical 
for low-mass disks in shearing potentials, but none are 
circular rings. They resemble the very winded, floccu- 
lent structures observed in real ETGs, and could trigger 
turbulent motions just like in th e star-forming spirals 
studied bv lElmegreen et all (|2003f) . But the Toomre pa- 
rameter remains above 1 (around 2) in these arms (Fig- 
ure [7|) , and the disk density remains below the tidal den- 
sity (see below and Figure [9|). No dense clump can then 
collapse and the star formation efficiency remains low. 
This is confirmed by the densest regions in these arms 
being short-lived and not collapsing into bound clumps, 

^ Note that particle noise and grid discretization make tlie veloc- 
ity dispersion uneasy to measure in a differentially rotating disk, 
our 15 km s~^ value is thus an upper limit, but careful tests have 
ensured that the actual dispersion is not lower than 10 km s~^. 



as shown by the evolution of the gas density over a short 
period (Figure [8]). 

The velocity dispersions of our modeled gas disk in 
a red ETG are thus likely driven by one or several re- 
alistic processes including non-axisymmetric sheared in- 
stabilities, which does not seem to be a numerical arti- 
fact. The resulting turbulent pressure can provide the 
support against star- forming axisymmetric instabilities. 
Observations of relatively high velocity dispersions in 
the non-star-forming ou ter HI disks of spiral galaxies 
(jTamburro et al.1 120091 ) . with level of velocity disper- 
sions close to that of star-forming disks, suggest that 
some processes do drive the turbulence^ without trigger- 
ing star formation, and could likely act in the internal 
regions of early-type galaxies to stabilize their gas disks 
in the morphological quenching mechanism. 

4.3. Star formation and resolution 

In the red ETG, the stability of the gas disk against ax- 
isymmetric perturbations prevents the formation of very 
dense gas clouds, for which high star formation efficien- 
cies would result. We indeed assumed in our simulations 
that the star formation rate depends on the local gas den- 
sity following a Schmidt law, which is realistic for various 
triggering mechanisms (pure gravita tional instabilities , 
turbulent distribution of the density (|Elmegreen Il2002f) . 
cloud-cloud collisions ([Smith I [1980( 1). and we impose a 
gas density threshold for star formation at space densi- 
ties of 0.03 Mqpc^'^, which is about one atom per cubic 
centimeter. This is a conservative assumption: for a gas 
disk of height 200 pc, it corresponds to a surface density 
of 6 Mqpc~^: this is the minimal density for cool diffuse 
atomic clouds formation, which is only the first step to- 
wards star formation. The critical density threshold for 
star formation is generally above this minimal density 
(see the review i n ,Elmegreen 2002. and observations by 
iKennicutt II1989D . 

We have performed a test simulation with a star forma- 
tion threshold set to 0.25 Mqpc"'^. The star formation 
history and color evolution are globally unchanged (Fig- 
ure |4|) , and the colors are even slightly redder between 
9 Gyr and 11 Gyr (with a maximum B — V — 0.90 at 
t — 11 Gyr). In the last part of the simulation, the galaxy 
experiences the same kind of bluing in both simulations. 
We also ran a simulation with no threshold for star for- 
mation, where the galaxy still becomes red during the 
quiescent phase (see Appendix A. 4 for more details). 

That actively star-forming instabilities do not arise 
during the red early- type phase is likely a robust result, 
not only because Q > 1, but also because the density 
of the gas disk remains below the limiting tidal density 
for the total galactic potential (Figure [9]). This means 
that even the densest regions of the gas disk, in the spiral 
armlets, cannot collapse because the tidal field exerted 
by the large-scale potential exceeds their internal grav- 
ity. The inefficient star-formation of the gas disk in the 
early-type galaxy is also consistent with this prescription, 
and the evolution of the gas density on short timescale 
(Figure [8|) shows that no bound clouds form even in the 
densest parts of the spiral armlets. 

* The velocity dispersion in these outer HI disks could however 
be mostly thermal (lElmegreen fc Parravano .,1994i) . 
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Fig. 9. — Ratio of gas density to the tidal density in the three phases of evolution, at t = 3.5, 9 and 13.7 Gyr respectively (for the 
simulation without the final major merger). The box side is 30 kpc. The early and late phases show extended regions where the ratio 
is above unity, so self-gravity wins over tidal stripping and allows the growth of perturbations. The MQ phase shows no such regions, 
explaining the absence of bound clumps in this phase. 



The ~ 250 pc scale-height of the stable gas disk in 
our simulation may be an artifact of the 130 pc spatial 
resolution. Had the disk been thinner, the gas would 
have been denser and less stable against star formation 
according to the tidal density criterion. However, such a 
very thin disk would require a much lower level of velocity 
dispersion, which would result only in small, low-mass 
star forming regions 

More generally, the mass resolution of a few IO^Mq 
makes it impossible to resolve or rule out the formation 
of dense, star-forming gas clouds with sizes of a few tens 
of parsec and masses of a few The formation 

of small clouds in stabilized gas disks needs to be in- 
vestigated in detail using simulations with a pc-like res- 
olution, which cannot yet be achieved in cosmological 
zoom-in simulations. Nevertheless, we can infer that any 
such star-forming instabilities would have been of small 
masses, sizes, turbulent speeds and star-formation rates 
compared to the classical star-forming regions of spiral 
galaxies. 

The MQ mechanism unveiled in our simulations re- 
mained unnoticed in previoi is studies of galaxy for ma- 
tion. This led for instance iKhalatyan et al.l (|2008f l to 
conclude that AGN feedback was likely necessary to ex- 
plain the reddening of ETGs. According to the analysis 
above, an accurate description of the color evolution of 
galaxies with gas disks relies on resolving star-forming 
clumps, which implies that the spatial resolution should 
be much better than one kpc, and the mass resolution 
should be at least ^ IO^Mq. The turbulent velocity dis- 
persion that acts to stabilize the disks should also be 
correctly simulated, which implies that shear instabili- 
ties such as spiral arms must be resolved. This requires 
the presence of a cold gas phase below lO^K. While such 
a phase is naturally present in our sticky particle models 
(see Appendix B), it requires a spatial resolution bet- 
ter than 100 pc in hydrodynamical simulations. A low 
spatial or mass resolution or a high temperature floor in 
the ISM modeling might prevent the formation of sub- 
structures and leave the gas disks smooth in all galax- 
ies, making it impossible to resolve dense star-forming 
clouds. 

5. DISCUSSION 

5.1. Morphological quenching in observed ETGs? 

A fundamental difference between the MQ mechanism 
and the classical quenching mechanisms like AGN feed- 
back and virial shocks is that red galaxies are not sup- 
posed to have their cold gas reservoirs removed or heated. 



A prediction if MQ plays an important role is that red 
early-type galaxies could contain significant amounts of 
cold gas, up to several percent of their baryonic mass. 
These cold gas reservoirs could sometimes reach masses 
and densities that would lead to efficient star formation 
and blue colors in spiral galaxies but do not form stars 
(or only inefficiently) when they are embedded in an el- 
liptical galaxy. 

In the Local Volume, the gas content of red ETGs 
(ellipticals and SOs) is often less than 1% of the bary- 
onic mass, if d etected at all (e . g.. IBregman et al. II1992I 
for atomic gas. ICombes et all 120071 for molecular gas). 
However, a large fraction of local ETGs are in the Virgo 
cluster or in dense groups, so they should have lost their 
cold atomic and molecular gas via ram pressure stripping 
(but see Tonnesen & Bryan 2009, for the effect of ram 
press ure on molecular c louds), harassment (Moore et al.l 
|1996[ ). strangulation (.Kawata fc Mulchaev, .2008.) . etc.. 
Recent studies suggest that HI de tections could actuall y 
be common among field ETGs (IMorganti et al\\200B). 
The gas fractions remain generally low (|Grossi et al. I 
120091 ) ■ and a significant fraction is often in low-density 
outer structures or debris from tidal tails rather than in- 
ternal disks (Oostcr loo et al I200a IMorganti et al.ll2006l : 
IOosterlo~t al.,.2007t) . but a variety of processes can be 
responsible for the low remaining gas masses (gas exhaus- 
tion, ram pressure and tidal stripping, strangulation) and 
the gas fraction could have been higher a few Gyrs ago. 

Nevertheless, dense gas is sometimes found in the in- 
ternal regions_ofETG^^OTenJn_t gas 
(e.g., ICrocker et all [2008I . [20091: lYoung et all [2009a b). 
Massive red ETGs can contain large amounts of dense 
gas in spite of their red colors, and only low levels of star 
formation are observed in these cases. A n oticeable ex- 
ample is the elliptical galaxy NGC 2320 (jYoung et al.l 
I2009bf) . This very red object {B - V ~ 1) contains 
about 5 X 10^ M0 of molecular gas in a disk of a few kpc 
of radius, not counting any possible atomic gas compo- 
nent, but its star formation rate remains at undetected 
levels, while such a high mass of molecular gas would 
turn any spiral galaxy blue. ESQ 381-47 is another red 
ETG that has an unusually low s tar formation e fficiency 
in its gas disk as re ported by iDonovan et al.l ([2009). 
IMorganti eFan (|2006f) show that NGC 3414 and 4278 
host large HI disks but only old stellar populations (see 
also Shapiro et al. in preparation). 

Although they do not dominate in the Local Universe, 
these examples of early-type galaxies with massive cold 
gas disks but red integrated colors and low star formation 
activities are likely candidates for the MQ mechanism. 
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If the morphological quenching mechanism plays a sta- 
tistically important role in forming the Red Sequence, 
a prediction is that red early-type galaxies would fre- 
quently host significant amounts of cold gas, of few per- 
cent of their baryonic mass. This would not necessarily 
be the case in the Local Universe where other environ- 
mental processes have exhausted the gas content. This 
could be more frequent at high [z > 1) or intermediate 
(z ~ 0.5) redshifts, when most ETGs are not in viri- 
alized clusters and dense groups yet, and the cold gas 
infall rates are higher. Recent st ellar population stud- 
ies of red galaxies at high-redshift (jMaraston et al.ll2CI09l : 
iTonini et al.l 120091 ) indicate that their colors are consis- 
tent with some residual star-formation activity. This 
could be compatible with the presence of gas disks form- 
ing stars inefhciently in a MQ mechanism, instead of 
complete r emova l of cold gas and termination of gas sup- 
ply. K riek et al.l (pOOi^) have recently obtained a near- 
infrared spectrum of one of the red ETG in their z ^ 2 
sample, which suggest a low but non-zero star formation 
rate. This could also be consistent with the presence of 
a morphologically-quenched gas disk. This red ETG has 
a stellar mass of 2x10^^ Mq and an estimated SFR of 1- 
3 M0 yr~^, which compares favorably to the MQ phase 
in our simulations, but at higher z. 

5.2. Observed blue ellipticals 

At the end of our simulation, we obtain a blue elliptical 
galaxy just after a major merger. This is also briefly the 
case at i ~ 6 Gyr. These are galaxies that have just been 
transformed into ellipticals by mergers and experienced 
a starburst at the same time. A galaxy that has under- 
gone a starburst can indeed remain blue for one or two 
Gyrs depending on the strength of the starburst and the 
photometric bands considered (e.g., Larson et al. 1980). 

Such blue elliptical galaxies are rare in the Local Uni- 
verse. Lee et al. (2008) find that 4% of ETGs at z < 0.1 
have blue colors, and most of these blue ETGs are low 
mass, dwarf ellipticals. There are nevertheless known 
cases of red ETGs with little on-going star formation 
but that have undergone a starburst in the last billion 
year: these E-l-A galaxies are rare at z = 0, proba- 
bly because galaxies that have undergone a starburst- 
triggering merger in the last couple of Gyrs are rare (e.g., 
iGenel et al] 120081 ). The blue elliptical at the end of our 
simulation results in fact from the combination of two 
different phenomena, a major merger at z — 0.15, and 
a high gas accretion rate between z — 0.4 and z = 
so that the final merger involves a gas-rich galaxy. Our 
simulation shows that this can happen in A-CDM cos- 
mology, although the rapidly decreasing major merger 
rate and gas accretion rate at z < 1 (e.g., [Kercs ct al. 
|2009() should make such cases infrequent, in agreement 
with the low fraction of E-l-A galaxies at low redshift. 
At higher redshift, the higher merger rate could increase 
the fraction of blue ellipticals and E-|-A galaxies, resem- 
bling for instance the blue elliptical at z ~ 1.5 following 
several minor mergers. 

The two blue ETGs in our simulations are short-lived: 
their recently formed stellar spheroid is quenching their 
star forming activity, and they are due to become red 
ellipticals in about 1 Gyr. When we continue the simu- 
lation after the z — instant, the final galaxy also turns 
back red. 



5.3. Revived blue disks? 

We found in our model a case of red spheroidal galaxy 
evolving into a blue galaxy with a massive disk. In the 
model without the final z ~ 0.2 major merger, this re- 
sults in a disk-l-bulge galaxy at z = with a grand-design 
spiral structure. However, our model does not support a 
"spiral rebuilding scenario" according to which the ma- 
jority spiral disks would be re-built after major mergers 
from tidal material fal ling back and/or fresh inf ailing gas 
(jHammer et al. ||2005() . because: 

(1) The reformation of a star- forming disk requires a 
high gas infall rate (~ 10 M© yr~^) during several Gyr, 
without any significant mergers during this period. This 
is above the average accretion rate expected from z ~ 1 
to z ~ for a typical ga laxy with the mass of the Milky- 
Way (jKeres et al.ll2009f). and above the t vpical observed 
gas infall rates, too (jSancisi et al. "1120081 ). Most galaxies 
would then be expected to accrete lower amounts of gas: 
in this case, the MQ mechanism remains active and the 
star formation is quenched, preventing the re-formation 
of a stellar disk and the return to blue colors. A "spiral 
rebuilding" mechanism below z ~ 1 then requires high 
rates of diffuse gas infall that seem unlikely to take place 
in most galaxies; the return of a red galaxy to blue colors 
is a possible mechanism but likely taking place only in a 
small fraction of galaxies. 

(2) The re-built disk galaxy keeps important signs of 
the high-redshift mergers, even in our model without the 
final major merger. A massive, red stellar halo, as large 
as ~20 kpc in radius, is present, amounting to about 
25% of the stellar mass at z = 0. The kpc-sized bulge 
amounts to 35% of the stellar mass, only 40% of the 
stellar mass hence being in the rotating disk. This disk 
galaxy rebuilt after a series of me rgers therefore differs 
from typical local spiral galaxies (jBinnev fc Tremaine I 
I1987I) . that have a significantly lower mass fraction in 
the rotating disk component. 

Furthermore, the transient blue elliptical phases in 
our models, either at z ~ 1.5 or z ~ 0, correspond to 
post-merger elliptical galaxies that are about to redden 
rapidly. Hence, neither of the blue ETGs found in our 
ACDM-based models have blue colors persisting more 
than 1 Gyr. They are not signs of sp iral galaxies be- 
ing re built after the mergers (but see iKannappan et al. I 
I2009D , and these blue elliptical phases, potentially E+A- 
like, result from the time needed for a stellar population 
to redden after the star formation is quenched. 

6. CONCLUSION 

We proposed a straightforward mechanism for the 
quenching of star formation in galactic disks once they 
become dominated by a stellar spheroid, which stabilizes 
the disk against fragmentation into bound, star-forming 
clumps. We demonstrated the natural occurrence of such 
morphological quenching in a cosmological simulation. 

For the purpose of a numerical study, we used a new ef- 
ficient technique for simulating the evolution of a galaxy 
in a cosmological context. The merger history and 
smooth accretion history are recorded for a given galaxy 
in a cosmological simulation, and the evolution of the tar- 
get galaxy is re-simulated at high resolution from z = 2 
to z = 0, including gas dynamics and star formation. 

We have studied the morphological and photometric 
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evolution of a field galaxy ending with a stellar mass of 
2x10^^ Mq at z = 0. This galaxy goes through three 
different phases of growth: (a) a blue phase from z ~ 2 
to 1, during which a stellar spheroid is built by a series 
of minor mergers, (b) a long red phase from z ~ 1 to 0.2, 
with smooth gas accretion at a high rate ~ 10 Mq yr~^, 
and (c) a blue phase after z ~ 0.2, following a major 
merger. Based on the analysis of the evolution of this 
galaxy, our main conclusions are as follows. 

Once a galaxy acquires a spheroid-dominated morphol- 
ogy, by a major merger, minor mergers or disk instabili- 
ties, it can spontaneously turn red. This occurs without 
gas removal or suppression of gas supply by other quench- 
ing mechanisms such as AGN feedback, virial shock heat- 
ing, gravitational quenching or ram pressure stripping. 

This morphological quenching in early-type galaxies re- 
flects the stabilization of the disk by the dominant pres- 
ence of a pressure-supported stellar spheroid, replacing 
the stellar disk. The stabilization of the gas disk is a 
result of two effects: the steep potential well induced by 
the spheroid and the removal of the stellar disk compo- 
nent that used to enhance the self-gravity of perturba- 
tions in the disk. Early-type galaxies can thus be red 
and dead even in the presence of a gas disk, either left 
over after mergers or accreted later from the cosmologi- 
cal reservoir. The MQ mechanism should be equally ef- 
ficient if massive stellar spheroids are formed by another 
process than mergers, for instance giant disk instabilities 
and clump coalescence at high redshift (Noguchi 1 19991: 
Elmegreen et al.l2008HDekel et al.ll2009bl : iCeverino et al l 
2009lV 

In the case of merging galaxies with a higher gas frac- 
tion than what we studied here, th e final disk to spheroi d 
ratio of the remnant is increased (|Hopkins et al.l [20091 ) . 
so that the effect of MQ would be weaker. However, 
most of these remnants are bulge-dominated, so that we 
can still expect the bulge to have a stabilizing effect, un- 
less the disk is particularly massive. Generally speaking, 
we propose that a galaxy with most of its stars in bulge 
and/or a thick disk might experience some morphologi- 
cal quenching, which could not only apply to ellipticals 
but also to SO galaxies. Future simulations will probably 
allow us to study this in more details. 

We also find that a red-and-dead galaxy can turn back 
into a blue star-forming galaxy. The gas disk can frag- 
ment to clumps and form stars again if it becomes dense 
enough to overcome the stabilizing effect of the spheroid. 
This requires a long phase of smooth gas accretion at a 
high rate, which is possible even at z < 1 but only for a 
small fraction of the galaxies. Such a galaxy that turns 



blue again is expected to have a spheroid of mass compa- 
rable to the disk mass, so it is not a likely progenitor of 
Milky Way-like spirals. Our models indicate that blue el- 
lipticals and E-fA galaxies are more likely to be recently 
formed spheroids that are about to turn red than red 
merger remnants that are about to re-create a massive 
blue disk. 

The natural process of morphological quenching pro- 
vides a simple explanation for the m ost distinct feature 
of th e observed galaxy bimodality (jPekel fc Birnboiml 
l2006f ). namely the fact that the color evolution of a 
galaxy is tightly linked to its morphological evolution. 
The MQ can help quenching the star formation in galax- 
ies on the Red Sequence. It should be naturally working 
in concert with other quenching mechanisms, depend- 
ing on the mass and environment of each galaxy. But 
a unique feature of MQ is the fact that it can explain 
the presence of red-and-dead galaxies in halos below the 
critical mass for virial shock heating, ^ IO^^Mq, where 
the other quenching mechanisms fail. 

Current observations of large atomic and molecular gas 
reservoirs in local early-type galaxies, with red colors and 
low star formation rates, suggest that the reddening of 
elliptical galaxies cannot be solely a result of gas removal 
or suppression of gas supply. Morphological quenching 
could thus contribute to the formation of the modern 
Red Sequence. The presence of quenched gas disks in red 
ETGs may be more frequent at higher redshift, where gas 
fractions and gas accretion rates are typically higher. 
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APPENDIX 

APPENDIX A. SIMULATION OF GALAXIES IN THEIR COSMOLOGICAL CONTEXT 

In the following subsections, we describe the initial cosmological simulation, the extraction of merger and accretion 
histories, and the technique used to perform the high-resolution re-simulation and to analyze it. 

A.l. Cosmological simulation 

The cosmological simulation has been performed with the Adaptive Mesh Refinement code RAMSES (|Tevssier I 
|2002| ). The simulation box contains 512^ dark matter particles and has a comoving length of 20 h~^ Mpc. This gives 
a mass resolution of 6.9x10^ M©, so that a Milky Way type halo is made of a few 10^ particles. The cosmology is set 
to A-CDM with fi™=0.3, nA=0.7, Ho=70 km s'^ Mpc'^ and (T8=0.9. 
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A. 1.1. Halo finding 

Th e dark matter halos are detected in the cosmological simulation using the HOP algorithm (lEisenstein fc Hutl 
Il998f ). with (5pcak=240, ^saddie=200 and (5outoi— 80. The minimal number of particles per halo is fixed to 10 (which 
corresponds to a minimal halo mass of 6.9x10^ M©). 

Particles that do not belong to a gravitationally bound halo will also be taken into account; we will refer to them as 
"diffuse" mass accretion when they are accreted by the studied halo. 

A. 1.2. Selection of the studied halo/galaxy 

We call "main halo" the halo (and galaxy) which we want to study and re-simulate, chosen in the final z = 
snapshot. When the studied halo is located in a rather empty environment (i.e. it does not belong to a cluster, and 
no more massive halo is found within 2 Mpc at z = 0), we do not take into account the large-scale tidal field when we 
re-simulate the history of the halo. This is the case for the halo studied in the present paper (see right panel of Figure 

A. 1.3. Merger and accretion histories extraction 

The method used to extract the merger and accretion histories consists in following the main halo (more exactly the 
most massive progenitor of the final halo) from z — 2 to z = and recording all the halos - the merger history - and 
diffuse particles - the accretion history - crossing a fixed spherical boundary around it. This boundary is a sphere of 
radius equal to the final virial radius of the main halo at 2; = 0, which is defined as: 



3A/t. 



4tt X 200pc 



where Aftot is the mass of the main halo at z = 0. 

At z = 2, the radius of this sphere is larger than the initial main halo: it encloses not only this main halo but also 
some satellite halos and diffuse particles. The re-simulation will begin at z = 2 with all the halos and diffuse particles 
initially enclosed in the boundary. To this aim, we record their mass, position, velocity, and spin for the satellites, all 
in the rest frame of the main halo. Positions are expressed in physical units instead of comoving, but with velocities 
corrected for the Hubble flow: a radial velocity Vr = H{z) x d (where d is the distance to the center of mass of the 
main halo) is added to the velocity of diffuse particles and to the velocity of the center of mass of the halos. 

We then check each subsequent output of the cosmological simulation (snapshots are 60 Myr apart), and search 
for particles that are inside the boundary but have not been found there earlier-on: these are particles that have 
crossed the boundary since the last snapshot. After having undergone the same modifications as described previously 
(i.e. coordinate transformation from comoving to physical and addition of the Hubble flow), their mass, position and 
velocity are written in a catalog (together with the redshift of their entrance in the boundary) . 

As far as halos are concerned, we record a halo if at least one particle belonging to it has crossed the boundary, 
and write in the catalog the position, velocity, mass and spin of this halo when its first particle(s) have crossed the 
boundary. The center of mass of the halo is thus still outside the boundary, particularly in the case of a large halo, 
but this choice has the advantage of being straightforward and accurate. According to its orbital parameters, the halo 
then might not totally enter the boundary and merge with the main halo. This is not a problem for the re-simulation, 
where it will not merge with the main halo either and will simply make a fly-by. 

The last point is that we have to handle the case of particles that enter a companion halo while this one has already 
been written in the catalog (i.e. the halo has grown since some of its particles have first crossed the boundary and 
this halo still has not merged with the main halo). In this case, the particles are handled as diffuse accretion. A halo 
can then be surrounded by a cloud of diffuse particles. This approximation concerns only the "satellites of a satellite" 
merging simultaneously, which concerns only a small fraction of the total infalling mass. 

To summarize, the final products of this analysis are the following: 

• an initial conditions file: a list of all the halos and particles that are already inside the boundary (the final virial 
radius) at z = 2 

• a catalog: a list of all the particles that cross the boundary between z = 2 and z = and of all the interacting 
companion halos. 

A. 2. High-resolution re- simulation 

We now replace each halo of the cosmological simulation by a galaxy made up of gas, star and dark matter particles, 
and each diffuse particle of the cosmological simulation with a blob of lower-mass, higher-resolution gas and dark 
matter particles. They start interacting with the main galaxy following the orbital and spin parameters given by the 
cosmological simulation. 

The code we use for the re-simulation is the Particle-Mesh code described in Bournaud & Combes (2002, 2003). 
The density is computed thanks to a Cloud-In-Cell interpolation and the Poisson equation is then solved using Fast 
Fourier Transform techniques. Time integration is made using a leapfrog algorithm with a time step of 1.5 Myr. The 
maximal spatial resolution (and softening length) is 130 pc. The mass resolution is 1.4x10^ M© for stellar particles 
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Fig. 10. — Initial distribution of stars (top panel) and gas (bottom panel) for the main galaxy, seen face-on and edge-on (each panel is 
40 kpc X 40 kpc in size). The small overdensities that can be seen belong to the disk, and are the result of gravitational instabilities due 
to the high initial gas fraction. 

initially present in the galaxies, 2.1x10^ M0 for gas particles and stellar particles formed during the simulation, and 
4.4x10^ M© for dark matter particles. 

Gas dynamics is modeled with a sticky particle scheme with f3r=0.8 and /3t=0.7:these are the factors by which the 
relative radial and tangential velocities of colliding gas p articles are reduced after the collision. Star formation is 
computed with a Schmidt-Kennicutt law (jKennicutt Ill998t ) : the star formation rate is proportional to the gas density 
to the exponent 1.5, above a fixed threshold of 0.03 Mopc"^, which is about one atom per cubic centimeter. The 
efficiency of star formation is chosen so that at z = 2 the main galaxy woul d consume its gas with a timescale of 
1.5 Gyr, consistent with the star-formation timescale of high-redshift galaxies (jPaddi et al. I [20071 ). 

Feedback from supernovae explosions is not included. Indeed, the Schmidt-Kennicutt law, which is an observed 
relation, should already account for this effect. What is more, simulations have shown that feedback only has a 
minor role in quenching star formation (Kravtsov 2003), except maybe in the case of very strong starbursts where 
strong winds could be triggered. We have also chosen not to include feedback from AGNs. AGN feedback could be 
implemented in our re-simulation technique just as in cosmological zoom simulations, but in the present paper we 
present a quenching mechanism that can take place independently of AGN activity. 

A. 2.1. Model galaxies 

Each halo recorded in the dark matter-only cosmological simulation is replaced with a realistic galaxy (Figure 
10), made of a disk (gas and stars), a stellar bulge and a dark matter halo ^. The di s ks have a Toomre profile , 
the bulge is a Plummer sphere and the halo follows a Burkert profile (jBurkert 1 119951 : iSalucci fc Burkerti l2000f) . 
This halo profile has been chosen because "real" galaxies most probably have halos with a cor e rather than with a 
cusp. The fact that dark matter halos have a fiat core is widely observed in the Loca l Universe ([Borriello fc SaluccTI 
I200H ISalucci II2001I: Ide Block et al. II2003I : lOh et al. 1120081 iTrachternach et al. Il2008f l and high redshift observations 
(jForster Schreiber et al.l 120061 ) suggest that the mass in the inner regions of galaxies at z^2 is essentially baryonic 
(which might not be the case if dark matter halos had a steep profile). 

The core radii of our dark matter halos are chosen according to the scaling relations proposed by Salucci & Burkert 
(2000) and the halos are truncated at their virial radius. 

The total mass of the galaxy is divided in 17% of baryons and 83% of dark matter (the dark matter mass is the mass 
of the halo found in the cosmological simulation) . The gas content in the disk varies according to the redshift at which 
the galaxy is introduced in the simulation and according to its tota l mass: for small galaxies (Mhaio < IO'^-'^Mq) the 
gas disk represents 30% of the stellar disk ([Roberts fc Havnes 1119941 ) . whatever the redshift, and for massive galaxies 
it is 30% at "high" redshift (z > 0.8) and 15% at "lo w" redshift (this is in order to mat ch the very high gas fractions 
observed at z ~ 2 for instance bv lDaddi et al. l[2008l and the low redshift values of iRoberts fc Havnes 1994). 

To decrease the computing time while keeping a realistic cosmological context, a different galaxy is not created for 
each halo of the cosmological simulation, but instead a set of model galaxies spanning the whole range of masses needed 
(and closely sampling the masses of the halos) is created (Table [2| shows the characteristics of these model galaxies) . 
When a halo has to be introduced in the simulation, the model galaxy that is the closest in mass is chosen among this 
set (note that the most massive companions are exactly reproduced). 

Each model galaxy is evolved in isolation for 500 Myr before being introduced in the simulation (star formation 
is shut down during this period so that the initial gas fraction is conserved). This allows the galaxies to acquire a 
realistic structure (e.g. spiral arms, bars...) before beginning to interact with one another, so that we can be sure 
that all the processes observed during the interactions are due to the interactions themselves and not to the growth of 
non-axisymmetric structures from axisymmetric initial conditions (which is not what we are interested in). 

^ Note that no hot gas halo is included 
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Fig. 11. — Large scale view of the gas distribution at the beginning of the simulation (the panel is 440 kpc x 440 kpc in size). The disk 
of the main galaxy is seen edge-on and is surrounded by satellite galaxies and diffuse gas as prescribed by the cosmological simulation at 



TABLE 2 
Characteristics of model galaxies. 





Gl 


G2 


G3 


G4 


G5 


G6 


G7 


G8 


G9 


GIO 




Dark matter halo 


Mass (Mq) 


1.7x10" 


5x10* 


1 X 10'^ 


2.5x10^ 


5x10^ 


7.5x10^ 


1x10"' 


2.5x10"' 


7.5x10"' 


4.1x10" 


^particles 


400000 


1136 


2272 


5681 


11363 


17045 


22727 


56818 


170454 


931818 


Rtrunc (kpc) 


38.63 


8.21 


10.34 


14.03 


17.68 


20.24 


22.27 


30.23 


43.60 


102.41 


Rcoro (kpc) 


25.80 


1.10 


1.40 


2.20 


3.00 


3.60 


4.10 


6.40 


10.80 


20.70 


Stellar disk 


Mass (Mq) 


2.5x1010 


SxlO'^ 


1.5x10** 


3.8x10** 


7.7x10** 


1.1x10^ 


1.5x10^ 


3.8x10^ 


9.18x10^ 


5.5x10"' 


Nparticlcs 


170078 


523 


1047 


2618 


5236 


7855 


10473 


26183 


62615 


376915 


I^.trunc (kpc) 


6.00 


0.32 


0.45 


0.71 


1.01 


1.24 


1.43 


2.26 


3.91 


9.15 


Rcarac (kpc) 


3.00 


0.16 


0.22 


0.36 


0.50 


0.62 


0.71 


1.13 


1.96 


4.57 


h (kpc) 


0.17 


0.01 


0.01 


0.02 


0.03 


0.03 


0.04 


0.06 


0.11 


0.26 


Gas disk 


Mass (Mq) 


7.5x10^ 


2.3x10^ 


4.6x10^ 


1.1x10** 


2.3x10** 


3.5x10* 


4.6x10** 


1.1x10^ 


2.8 X 10^ 


8.3x10'-' 


^particles 


357165 


1099 


2199 


5498 


10997 


16495 


21994 


54985 


131493 


395760 


Rtrunc (kpc) 


9.85 


0.52 


0.74 


1.17 


1.66 


2.03 


2.35 


3.71 


6.43 


15.03 


Rcarac (kpc) 


8.13 


0.43 


0.61 


0.97 


1.37 


1.68 


1.94 


3.06 


5.31 


12.42 


h (kpc) 


0.07 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.03 


0.05 


0.12 


gas fraction 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.30 


0.15 


Stellar bulge 


Mass (Mq) 


2.6x10'' 











1.4x10*^ 


2.1x10^ 


2.8x10'' 


7.3xl0« 


2.9x10" 


1.7x10^" 


^particles 


18897 











10 


15 


20 


52 


20871 


125638 


Rtrunc (kpc) 


1.97 











0.33 


0.40 


0.47 


0.74 


1.28 


3.00 


Rcarac (kpc) 


1.11 











0.19 


0.23 


0.26 


0.42 


0.73 


1.70 



A. 2. 2. Diffuse mass accretion 

The catalog we have extracted from the cosmological simulation also takes into account diffuse accretion of dark 
matter particles onto the main halo. In the re-simulation, we replace each of these diffuse particles by a blob of gas 
particles (17% in mass) and dark matter particles (83% in mass). Each blob is thus made of 65 gas particles and 15 
dark matter particles. 

These particles are randomly distributed in a sphere. The radius of this sphere is chosen such that the mean density 
of the blob is lOOpc, which is chosen to be lower than the local density of the dark matter halo of the main galaxy so 
that the blob can be tidally disrupted when approaching the main galaxy: this allows the diffuse gas density field to 
be smooth, at least in its the densest regions (Figure [TT]). To respect this density criterion, the blobs introduced in 
the simulation before z = 1 have a radius of 2.5 kpc, the blobs introduced between z = 1 and z = 0.5 have a radius of 
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3.3 kpc and the blobs introduced after z — 0.5 have a radius of 5 kpc. 

The particles belonging to the blob are given no velocity dispersion since we expect the blob to be disrupted by the 
tidal field of the main galaxy anyway. 

A. 2. 3. Carrying out the simulation 

We start the simulation with the main galaxy (galaxy Gl in Table [2]) at z = 2 and place around it the satellite 
galaxies and gas blobs as prescribed by the initial conditions file (see Figure [TT] the result of this operation) . To place 
a satellite galaxy, we choose among the model galaxies the closest in mass, and put it at the right location with the 
right velocity. Its orientation is set assuming that the direction of the angular momentum of the disk is the same as 
the direction of the angular momentum of the dark matter halo detected in the cosmological simulation. 

Then, we read the catalog line per line and add in the simulation the other galaxies and the other blobs at the exact 
moment when they appeared in the cosmological simulation, with the exact position, velocity and spin extracted from 
the cosmological simulation. 

A. 3. Analysis 
A.3.1. Global mass evolution 

To study the evolution of the total mass of stars and gas in the main galaxy, we take into account all the particles 
located within a sphere of radius 25 kpc around the center of the galaxy. In the same way, to study the star formation 
rate of this galaxy, we only count stars that form within this sphere (i.e. we do not take into account stars forming in 
approaching companions that have not merged with the main galaxy yet). 

A. 3. 2. Stellar ages and colors 

Our simulation contains two types of stellar particles: those that were formed during the simulation (and for which 
determining an age is straightforward) and those that were initially present in the model galaxies we introduced in the 
simulation. We need to give an age to these stars in order to be able to compute galaxy colors, and for that we have 
used and tested different assumptions. Their common characteristic is that they assume that the first stars began to 
form 500 Myr after the Big Bang (that is around z = 10). The subsequent star formation histories have been tested 
(they concern not only the main galaxy, but also the infalling satellites for which we call Tin the moment when they 
are introduced in the simulation): 

• an exponentially decreasing SFR with a timescale Tgf given by T^/ — Tm — 500 Myr 

• an exponentially decreasing SFR with a timescale Tsf — 1 Gyr 

• a constant SFR between 500 Myr and 

• half of the stars formed at i = 500 Myr and the other half formed according to a constant SFR between 500 Myr 
and Tin 

Knowing the age of a stellar particle, w e use the spect ral evolution model P egase.2 to compute stellar color s 
assuming a Salpeter initial mass function (jSalpeter I [ToSSh from 0.1 to 120 Mq (|F ioc fc Rocca-Volmerange lll999( ). 
Each stellar particle is thus given a luminosity in B, V and K bands (Johnsson filters), all magnitudes being then 
computed in the Vega magnitude system. 

We assume a uniform solar metallicity in all our computations. Indeed, massive g alaxies at z y 2 are known to have 
a rela tively high metallicity already, even if somewhat sub-solar in their outer disks (|Erb et al. II2006I : iBournaud et al.l 
|2008() . and the B — V col or (on which we are going to focus) is not significantly sensitive to the metallicity (contrary 
to some other bands, e.g. iLarson et al. Ill980h . 

The integrated galaxy luminosity in various bands is then computed by summing up the stellar population spectra 
(we then consider all stars included in a sphere of radius 25 kpc around the center of this galaxy). Extinction is not 
taken into account, since the evolution of galaxies between the blue and red populations is known not to be caused 
by dust absorption effects; only some specif ic cases have their integrated color significantly reddened by extinction 
(|Alam fc Rvden II2002I : lUnterborn fc Rvden | [2008V Neglecting dust reddening is a conservative assumption regarding 
our main conclusion on the spontaneous reddening of elliptical galaxies. 

A. 4- Sensitivity to star formation threshold 

In the fiducial simulation shown in this paper, star formation is modeled with a Schmidt-Kennicutt law (|Kennicutt I 
Il998| ) above a threshold of 0.03 Mqpc^^. This is a conservative assumption since this threshold corresponds to the 
minimal density for cool diffuse atomic clouds formation (see also Section 4.3). To test the influence of this choice on 
our results, the same simulation has been run with a higher threshold (0.25 M0pc~^, which is possibly more realistic) 
and with a threshold set to 0, which is unrealistic but useful for comparisons. 

We find an overall similar evolution of the star formation rate in all cases (see left panel of Figure [T^ . Most of the 
differences arise during the morphological quenching phase, where, when no threshold is used, the star formation rate 
is twice higher with respect to what was found for the fiducial simulation. Similarly, setting a higher threshold reduces 
the star formation rate during this phase. As a consequence, the main galaxy experiences a reddening during the 
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Fig. 12. — Star formation rate (left panel) and corresponding color evolution (right panel) for three different star formation recipes : 
without any threshold, with a threshold set to 0.03 Mqpc~^ (our standard conservative assumption) and a threshold set to 0.25 Mqpc~^ 
(a possibly more realistic assumption). During the morphological quenching phase, all models exhibit a lower star formation rate and a 
redder color, even if the effect is stronger for higher thresholds. 

morphological quenching phase for the three different models, but this reddening is stronger for higher star formation 
thresholds (see right panel of Figure [T^ . 

Our results on morphological quenching (that were obtained using a very conservative assumption on the value of 
the star formation threshold) are thus robust. Using a higher (possibly more realistic) threshold would only increase 
the MQ effect. 

A. 5. Sensitivity to initial conditions 

To test whether or not our results depended on the initial conditions chosen (i.e. for instance gas fraction or Hubble 
type of the main galaxy Gl), we have run the exact same simulation replacing Gl with an elliptical galaxy. This 
galaxy is modelled as a Plummer sphere of stars (with a characteristic radius equal to 2 kpc and truncated at 6 kpc), 
with the same total baryonic mass as Gl (the dark matter halo is the same as for Gl). This allows us to test the 
influence of the gas fraction and of the Hubble type at the same time. 

In the first 1.5 Gyr of the simulation, we find a significant difference between this new simulation and the one 
presented so far. Indeed, since the initial galaxy now does not contain gas at the beginning of the simulation (Figure 
13), its star formation rate is much lower and thus its color is much redder (see Figure [HI). 

Yet, the galaxy acquires gas very quickly (mainly thanks to accretion from the filaments) and after 1.5 Gyr the gas 
fraction becomes remarkably similar in both simulations (see Figure [T5|) . Following this, the SFR evolution is the 
same, and thus the color evolution is also the same (see Figure [14]). 

In fact, the change in the initial conditions only affects the very beginning of the simulation both in terms of 
morphology and color. Our results are robust versus changes in the initial characteristics of the main galaxy since our 
mechanism to produce a red galaxy only requires: 

• a stellar spheroid, that is formed independently of the initial conditions (a disk will be destroyed by the first 
series of mergers, an elliptical will just stay elliptical and grow) 

• a gas disk that is here fueled by the filament and does not depend on the initial conditions either. 

APPENDIX B. THE STICKY PARTICLE SCHEME FOR GAS DISKS DYNAMICS 

Most of the gas mass in galac tic disks is in a cold atomic phase and molecular clouds rather than in warm and hot 
diffuse phases (e.g., [Brinks' 1990). T his rnedium is very h eterogeneous and its turbulent pressure dominates the thermal 
pressure (see reviews in Elmegreen 11200^ lBurkert1l2006[ ). It cannot t hen be describ ed as an ideal, continuous gas, and 
a description based on separate, interacting clouds is more adequate ([Larson II1981D . A particle-based hydrodynamical 
code, like SPH codes, could offer such a description only if dense and cold small-scale substructures are resolved 
together with their turbulent motions. This cannot be achieved with the 130 pc softening length of our simulations^. 
A temperature floor of a few 10'^ or lO** K would have to be imposed, as would also be the case for grid hydrodynamic 
codes (AMR) if a spatial resolution of a few parsecs cannot be achieved. The ISM modelled with such codes would 
then remain dominated by a high thermal pressure with a minimal sound speed of ^ 10 km s^^, and homogeneous on 
scales of at least 1 kpc. Such a description would be inadequate for star-forming gas disks, and would be even more 
problematic regarding the stabilization of a gas disk in an ETG, as the minimal sound speed imposed would artificially 
stabilize such disks independently of their real turbulent support. 

^ at fixed CPU time, an SPH simulation would actually have a lower spatial resolution and number of particles than our sticky particle 
simulations. 
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Fig. 13. — Evolution of the dense (p > 0.01 Mqpc ^) gas fraction (with respect to total baryonic mass) in a sphere of radius 25 kpc 
around the main galaxy, in the case where this main galaxy is a spiral at z = 2 and in a case where it is an elliptical. 




Fig. 14. — Color evolution of the main galaxy in the case where this main galaxy is a spiral at z = 2 and in a case where it is an elliptical: 
after 1.5 billion years, the color evolution is independent of the initial conditions. 

This motivates the use, in our Ral axy-scaled. re -simulations, of the sticky particle (SP) scheme, which is a realisti c 
description for a coUisional medium (IBrahidll977l ) or a cloudy, heterogeneous gas medium (Levinson & Roberts TOST). 
In our SP code, gas particles interact with each other through inelastic particle-particle collisions. The velocity 
dispersion of particles represents the turbulent speed of the modelled gas medium - not its thermal sound speed. 
The cross-section given to gas particles ensures about one collision per particle and per vertical crossing time and 
corre sponds to a filling factor in spiral disks of about 5%, which can be realistic for the ISM (|de Avillez fc Breitschwerd^ 
120041 ). Colliding particles bounce back, which models the turbulent pressure acting on gas substructures. These 
collisions are inelastic. We use a restitution factor /? = 0.75 for velocities (about 50% for the energy), because 
two gas clouds with similar masses and densities, surrounded by vacuum, that collide with each other would keep 
half of their initial kinetic energy in a kinetic form. The dissipated kinetic energy is converted into internal energy 
for each gas cloud (particle) and radiated away. With the collision rate of particles, this ensure that the velocity 
dispersion dissipates in about a vertical crossing time (in the absence of triggering source) , as found in high-resolution 
hydrodynamical sim ulations of turbulence (Mac Low 1999). The complete equations of the SP code can be found in 
iBournaud fc Combes I (|2002) and references therein. 

The SP scheme is not intended at being adequate on large cosmological scales. However, it could also describe 
the hot gas phase if collisions between particles are considered as the analogs of collisions between atoms: SP could 
then also reproduce thermal pressure. I t is neve rtheless most frequently employed to model the dynamics of the 
cool and cold ISM in galactic disks fe.g.. ISchwarz 1981; Noguchi 1999). An extensive study of cloud formation and 
disruption in SP simulat ions has been perform ed by iCombes fc GerinI (119851). Comparison of SP and SPH for disk 
dynamics can be found in lBournaud et al] (|2005l ). Parametric studies in lBournaud fc Combes I (|2002l appendix A) and 
IBournaud et al.l (|2005f l show that over a certain range, variations of the dissipation parameter f3 have limited effects on 
the global evolution of modeled gas disks as a stronger dissipation is compensated by a stronger gravitational heating. 
Motivation s for using an SP scheme ra ther than SPH at the scale of galaxies, in particle-based simulations, are further 
reviewed in lSemelin fc CombesI (|2002D . who show that the formation and evolution of disk galaxies in the cosmological 
context can be realistically studied with an SP code. 
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Fig. 15. — U-B and U-V color evolution. The shaded areas along the curves represent the range of colors obtained for different star forma- 
tion histories (Appendix . Both colors evolve similarly with time and in both cases a clear reddening is found during the morphological 
quenching phase. 

APPENDIX C. U - B AND U-V COLOR EVOLUTION 

We have studied the U ~B and U — V color evolution (see Figure fTS)) . We find results that are in good agreement with 
what we deduced from the B — V color evolution, in particular the reddening during the morphological quenching phase 
is a robust fact. For instance, U — B grows from —0.7 during the initial starburst phase to during the morphological 
quenching phase. 

U—B = is not a very red color for a galaxy (for instance. IWeiner et al. |[2005l define the red sequence as U—B > 0.1), 
but taking extinction into account could slightly increase this value and make it compatible with observations of red 
sequence galaxies, in particuler at high redshift. 
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